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A unified view on macroscopic thermodynamics and quantum transport is presented. Thermo- 
dynamic processes with an exchange of energy between two systems necessarily involve the flow 
of other balanceable quantities. These flows are first analyzed using a simple drift-diffusion model, 
which includes the thermoelectric effects, and connects the various transport coefficients to certain 
thermodynamic susceptibilities and a diffusion coefficient. In the second part of the paper the con- 
nection between macroscopic thermodynamics and quantum statistics is discussed. It is proposed 
to employ not particles, but elementary Fermi- or Bose-systems as the elementary building blocks 
of ideal quantum gases. In this way, the transport not only of particles, but also of entropy can 
be derived in a concise way, and is illustrated both for ballistic quantum wires, and for diffusive 
conductors. In particular, the quantum interference of entropy flow is in close correspondence to 
that of electric current. 
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I. INTRODUCTION 

Usually, macroscopic thermodynamics is derived on 
the basis of statistical mechanics and hence considered to 
be of essentially mechanical origin. This review provides 
a complementary point of view. It takes thermodynam- 
ics as the starting point and presents an elementary and 
transparent description of the quantum transport of par- 
ticles and entropy, both in the ballistic and the diffusive 
limit. Many results, scattered throughout the literature, 
are re-derived within a single, coherent approach. 

As a heritage of the classical view on physics, parti- 
cles interacting via certain force fields are considered as 
the elementary building blocks of matter. The evolution 
of quantum mechanics has forced us to accept that this 
statement can only be kept valid, if the meaning of the 
word 'particle' is considerably modified: rather than be- 
ing a point-like classical object with well defined position 
an momentum, particles of any kind have been under- 
stood as the quantized excitations of different quantum 
fields. In this way a unified description of fields and mat- 
ter has been achieved, which is known as 2 nd quantiza- 
tion. On the other hand, physics education still has to 
start out in the classical world. This leaves a deep imprint 
in our mind and language. When we want to proceed into 
the regime of quantum physics this imprint is not easy to 
overcome. This becomes particularly obvious in the field 
of quantum transport. 

The transport of physical quantities like energy, mo- 
mentum, angular momentum, electrical charge or en- 
tropy (E, P, L, Q, or S) can be described in terms of the 
transport of suitably chosen particles or quasiparticles, 
which carry the other physical quantities in the sense 
that with each particle a certain amount of E, P, L, Q, 
or S is associated, which is transported along with the 
particles. This means that the transport of particles is 
unavoidably connected to the transport of other physi- 
cal quantities. Importantly, only those physical quanti- 
ties can be transported for which a balance is possible, 



which tells us which amount of this quantity has left at 
time t\ a certain volume element V\ in space and has ar- 
rived at a time ti in another volume element V2 in space. 
The prototype of such balances are those of amounts of 
a substance, of goods or of money. It is suggested to 
call the physical quantities allowing similar operations 
balanceable or substance-like. A local density x, and a 
local current density jx can be associated with each bal- 
anceable quantity X. This property allow to formulate 
conservation laws in term of a continuity equation for 
each balanceable quantity (see below). 

Note that balanceable quantities are not necessarily 
conserved: important examples are entropy, or the spin. 
The total angular momentum is of course conserved, but 
the spin of the moving particles under consideration can 
be transferred to other systems, i.e., by spin-dependent 
scattering processes. Entropy can be generated, without 
extracting it from another system. Particles like pho- 
tons, phonons or excitons can be created or annihilated, 
provided that the canonic conservation laws for energy, 
momentum, and angular momentum are obeyed. 

Although the description of transport in terms of the 
motion of particles can be very easily visualized, it is hard 
to avoid the traps of classical physics in doing so. The 
reason is that quantum properties, in particular the in- 
distinguishability of identical particles, have no counter- 
part in the classical world, although they clearly show up 
in the thermodynamic and transport properties of mat- 
ter at the macroscopic level. The goal of this article is 
to formulate a description of transport processes in sys- 
tems of indistinguishable particles, which does not con- 
tain elements incompatible with the statistical concepts 
of quantum physics. 

The building blocks of this description are the elemen- 
tary Bose- or Fermi-Systems introduced below. Inter- 
estingly, the nature of these systems is microscopic and 
macroscopic at the same time. They are not point-like, 
but (in the same sense as a wave function) spatially ex- 
tended, they obey the laws of macroscopic or phenomeno- 
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logical thermodynamics, and the quantum mechanical 
average value of their particle number can be small - 
i.e., down to one or even less. In conventional language 
these systems are termed 'single particle states' and can 
be populated by one (in the case of fermions) or many 
(in the case of bosons) particles. The advantage of intro- 
ducing such systems - which have no place in a classical 
world - is the following: they allow to avoid a terminol- 
ogy, according to which non- or weakly interacting quan- 
tum particles occupy the states of a 'one-particle system', 
which is described by the conventional Schrodinger equa- 
tion. Although this terminology is well established, it 
leads to a confusion of the concepts of 'state' and 'sys- 
tem', which is very disadvantageous if one tries to formu- 
late a consistent thermodynamic description of systems 
with identical particles. 

In the course of the development of thermodynamics 
it became apparent that the traditional concept of heat 
had to be split into two more abstract concepts: 1 one, 
entropy, which is specific for thermal phenomena, and 
a second one, energy, which is relevant in all branches 
of physics. The terms 'heat' and 'heat current' survived 
in modern physics as a concept characteristic for pro- 
cesses, namely the amount of energy transferred from 
one system into another together with a given amount of 
entropy. This conception of heat, however, is problem- 
atic: despite the fact that heat cannot be balanced, 2 nor 
linked to quantum states, 3 it often competes with the 
much more powerful concept of entropy. Even today, 4 
entropy is less popular than heat, possibly because of 
persistent tradition, and its incompatibility with classi- 
cal mechanics. On the other hand, entropy can be han- 
dled very easily at the macroscopic level, as it behaves 
as an analogue of electric charge (except for the property 
of conservation). In this article, we take the existence of 
entropy as a starting point for thermodynamics - in the 
same way, as one takes the existence of electric charge as 
starting point for electricity. Like in electrodynamics, we 
don't ask what electric charge or entropy actually 'are', 
but take them as fundamental concepts, which prove use- 
ful in the quantitative description of electric or thermal 
phenomena, respectively. An intuition for both charge 
and entropy can be developed only via the many exam- 
ples, where we see them 'at work'. Our approach allows 
a natural integration of the ideas of quantum physics. 

This article is organized as follows: in section II we 
first formulate thermodynamics in a self-contained way 
that is appropriate for the investigation of transport pro- 
cesses. In section III we present the simplest of all trans- 
port theories, the drift-diffusion model, which has the 
great advantage of providing a simple intuitive picture 
of diffusive transport. In this approach particles and en- 
tropy are treated on the same footing. In section IV, we 
derive quantum statistics by combining thermodynam- 
ics with elementary ideas of quantum physics. In sec- 
tion V we introduce elementary Fermi-and Bose-systems 
as the elementary building blocks of quantum gases, and 
derive their thermodynamic equations-of-state. In sec- 



tions VI and VII we apply these equations-of-state to 
ballistic transport in one-dimensional quantum wires in 
the Landauer-Biittiker approach. In section VIII we show 
that the very same equations-of-state can be applied to 
generalize the drift-diffusion model in a way, which is 
equivalent to the Boltzmann equation in relaxation time 
approximation. From this perspective, there is no funda- 
mental difference between classical and quantum trans- 
port. In section IX some implications of our approach 
are discussed - in particular it is shown that also quan- 
tum interference can be included into the discussion of 
thermal transport phenomena. 



II. THERMODYNAMICS 

As pioneered by Massieu and Gibbs, thermodynamics 
can be based on the postulate that the static properties 
of any physical system with r independent variables can 
be compressed into certain functions of these variables, 
which are called thermodynamic potentials. The most 
familiar thermodynamic potential is the energy E when 
expressed as a function of the independent extensive vari- 
ables of the system. In the case of simple fluid or gas, 
the independent extensive variables are the entropy S, 
the volume V and the particle number TV. Assigning val- 
ues for a set of independent variables, e.g., for {S, V, N}, 
specifies a certain state of the system. The total differ- 
ential of the function E(S, V, N) can then be written in 
the form 



dE = TdS -pdV + ndN , 
where the absolute temperature T is defined as 



T(S,V,N) : 
the pressure p as 

- P (S,V,N) 



dE{S, V, N) 
dS 



dE(S, V, N) 



dV 



and the chemical potential \x as 



H(S,V,N) = 



dE{S,V,N) 
dN 



(1) 



(2) 



(3) 



(4) 



Equations 2, 3, and 4 can be called the caloric, ther- 
mal and chemical equation of state (EoS) of the system, 
and they define the to S, V, and N thermodynamically 
conjugate variables T, p, and fi. The EoS constitute a 
complete characterization of a specific system, and their 
knowledge is (up to an integration constant) equivalent 
to that of E(S,V, N). In mechanics or electrostatics a 
potential serves the purpose of combining three force or 
electric field components into a single function (the po- 
tential energy). In the very same way the energy (and 
also other thermodynamic potentials) combines the infor- 
mation contained in the three EoS into the single function 
E(S,V,N). 
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Equation 1 is called the Gibbs fundamental form 
(GFF) and has a very simple physical meaning: the in- 
tensive quantities T, —p, and [i tell how much energy 
has to be added to or removed from the system, if the 
extensive quantities S, V , and N are changed. In many 
cases a change of the extensive variables corresponds to a 
transport process: if the system is heated by connecting 
it to an external reservoir for energy and entropy while 
V and N are kept constant, the amount AE = TAS 
of energy together with the amount AS of entropy has 
to be transported from the reservoir into the system. If 
particles are added to the system at constant S and V 
from an external container, in addition to the amount 
AiV of particles the amount /J.AN of energy has to be 
transferred from the reservoir to the system. 

Depending on the specific problem under considera- 
tion, it is often convenient to use a different set of in- 
dependent variables by exchanging any of the extensive 
variables {S,V,N} with its thermodynamically conju- 
gate partner. If, e.g., {T, V, N } are chosen as set of in- 
dependent variables the corresponding thermodynamic 
potential is the free energy F(T, V, N) = E(T, V, N) - 
T ■ S(T,V, N). In the following the grand canonical, or 
Landau potential 

K(T, V, (i) = E — TS — fiN (5) 

will play a primary role, which is adapted to {T, V, /i} as 
set of independent variables. 

The fact that the EoS can be derived from a thermo- 
dynamic potential implies that they are not independent 
of each other. As the mixed 2 nd -order partial derivatives 
of continuously differentiable functions are equal, there 
must exist many relations between the partial derivatives 
of the different EoS. With {T, V, N} as set of independent 
variables one finds, e.g.: 

dS(T, V, N) _ _ d 2 F(T, V, N) 
dV ~ dVdT 

d 2 F(T,V,N) = dp(T,V,N) 

dTdV dT ' U 

Relations of this type are called Maxwell- Relations. 9 

The final feature of thermodynamics needed for 
the present work is the homogeneity postulate: 10 
it is assumed that the thermodynamic potentials 
Y(Xi, . . -Xj, Cj+ii ■ ■ • j £r) °f an y system have to be ho- 
mogeneous in the extensive variables {X±, . . . , Xj}, i.e. 

A Y(X\ , . . .Xj . . . . , £ r ) 

= Y(\Xi,...,\Xj,£j + i,...,£ r ), (7) 

where ■ ■ • , £r} are the independent intensive quan- 

tities, and A is an arbitrary dimensionless scaling factor. 
Equation 7 implies the Euler or homogeneity relation 

E = TS-pV + pN. (8) 

The physical meaning of the homogeneity 11 is quite fun- 
damental: it implies that all physical systems obey a scal- 
ing relation, which expresses the relations between their 



state variables in a way that is independent of the 'size' 
of the system, i.e., all properties of the system can be 
expressed by relations between the intensive quantities 
e, s, ?i, T, and fi . 

The homogeneity of E(S,V, N), for example, means 
that E can be written as 

E(S,V,N) = V-e(s,n) , (9) 

where e = E/V is the energy density, s = S/V the en- 
tropy density, and n = N/V the particle density. This 
implies that the function e(s, n) represents a reduced ther- 
modynamic potential, which still contains all thermo- 
dynamic information about the system, except its vol- 
ume. The corresponding reduced Gibbs fundamental 
form reads 

de = T ds + \xdn . (10) 

Equally well one can rewrite Eq. 9 as 

E(S,V,N) = N-e(s,v) , (11) 

where e = E/N is the energy per particle, s = S/N 
the entropy per particle, and v = l/n = V/N the vol- 
ume per particle. Then the function e(s,v) represents 
another reduced thermodynamic potential, which con- 
tains all thermodynamic information about the system, 
except its particle number. The corresponding reduced 
GFF reads 

de — T ds — pdij , (12) 

which is the form implicitly considered in many text- 
books, after stating that N is assumed to be fixed. The 
latter choice is preferred in physical chemistry, where it 
allows to express the specific properties of a substance in 
a given aggregation state in a way that is independent of 
the amount N of the substance. 

If we choose e(s,n) as reduced thermodynamic poten- 
tial, we can still apply the formalism of Legendre trans- 
forms to exchange the independent variables, i.e., s with 
T and n with fi. If we do so, we find using Eq. 7 

-p(T,n) = e(T, fi) - T • s(T, fi) - fj, • n(T, fi) (13) 

as the corresponding reduced thermodynamic potential. 
The differential of p(T, /j) represents a reduced fundamen- 
tal form, which is also known as the Gibbs-Duhem rela- 
tion 

dp = s dT + nd/j , (14) 
with the two equations of state 

8{T,fi) = dT and n(T,n) = — — — . (15) 

Using Eq. 8, we see that the reduced thermodynamic 
potential —p(T, ji) is equivalent to the Landau potential 

K(T,V,ii) = -Vp{T,ii) . (16) 



FIG. 2. An entropy current flowing through two isothermal 
surface elements with temperatures T\ and T2 < T\. The 
total rate of entropy production between the two surfaces be- 
comes negligible compared to the entropy currents through 
the surface elements, as the distance d and the temperature 
difference Ti — T2 go to zero. 



FIG. 1. A container with a gas of (quasi)-particles can be 
decomposed into small subvolumina with arbitrary size. Each 
subvolume represents another realization of the system 'gas', 
which continuously exchanges energy, entropy and particles 
with its neighbors. In presence of a gradient of T or jj, each 
subvolume can be still considered to be in local equilibrium, 
provided that its volume is not smaller than A 3 (A is the mean 
free path between scattering events, see Sec. III). The state 
of each subvolume is then characterized by the local values of 
T and /x. 

The choice of {T, V, /j,} as independent variables per- 
fectly matches the needs for a description of an impor- 
tant class of non-equilibrium situations. If we decompose 
a macroscopic solid, liquid or gas into small volume ele- 
ments as it is illustrated in Fig. 1, assuming that each of 
these volume elements is in local thermodynamic equi- 
librium, we can try to model situations where the local 
temperature and/or the local chemical potential are spa- 
tially varying. In this case the gradients of T and /1 are 
known to lead to the diffusion of energy, entropy and 
particles. 

It is appealing that also external helds like the grav- 
itational field or the electrostatic field can be built into 
the local thermodynamics. In a system with electrically 
charged particles we have one more extensive variable, 
i.e., the electric charge Q, which provides an extra term 
in the Gibbs fundamental form (see Eq. 1): 

dE = TdS -pdV + iidN + <t>dQ . (17) 

The to Q thermodynamically conjugate variable is the 
electrostatic potential <j>, which determines the electro- 
static contribution to the energy required for a local in- 
crease of the charge density. 

Usually charge and particle number are connected by a 
characteristic constant of the system, i.e., the charge per 
particle q. 12 In these cases charge and particle number are 
not independent, but proportional: Q = qN, implying 
that we can combine the last two terms in Eq. 17 

dE = TdS ~pdV + fidN (18) 



where 

P- ■= M + # (19) 

defines the electrochemical potential. For charged parti- 
cles it is p, and not /i, which quantifies the energy changes 
required for adding or removing particles. Hence p, and 
not just fi enters all thermodynamic relations for systems 
of charged particles. 13 

So far we have not exploited any of the conservation 
laws. 14 In general, a balanceable quantity X has to obey 
a continuity equation 

9X( ^ +V-j x (t, r) = Tix (t, r) , (20) 

where x is the local X-density, jx the AT-current den- 
sity and Yix the A"-generation rate per volume. For con- 
served quantities vanishes. Applying the continuity 
equation for E, 5, and N we can rewrite the reduced 
GFF (Eq. 10) in a very intuitive form: the rate of energy 
transfer e^lt^^nit^rj) from one side of an infinitesi- 
mally small surface element in space to the other side 
reads 

de(s(t,r),n(t,r f )) ds _ dn 

dt = T dt + »df (21) 

For a stationary flow, where all current densities are con- 
stant in time, the time derivatives of e, s, and n can be 
replaced by the corresponding current densities 

3 b = Tj s + p,j N , (22) 

because the total generation (or annihilation) rate of en- 
tropy and particles on the surface separating two volume 
elements vanishes, as its volume is zero. This is illus- 
trated in Fig. 2. 

Equation 22 provides a general relation between the 
energy current density and the current densities of the 
other independent balanceable quantities of the system. 
It is as fundamental as Eq. 1, and should hold for any sys- 
tem in local thermodynamic equilibrium, for which the 
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average velocity (v) is negligible. If local equilibrium is 
maintained in systems with non-negligible average veloc- 
ity, the (thermodynamically) conjugate variable pair ve- 
locity (v) and momentum (P) gives rise to another term 
in Eqs. 1 and 22 - including this term would lead to hy- 
drodynamics, and is beyond the scope of this article. For 
the diffusive transport to be studied in the next section 
(v) is usually so small that this additional term can be 
neglected. 

On a surface of constant T and p, an analogous expres- 
sion holds for the currents through this surface: 



fiI N 



This expression tells the strength of the energy current 
that is 'carried' by the currents Is and In, namely T Is, 
and ft, In, respectively. 15 The term T Is in this relation 
is usually called the 'heat current', but we have to stress 
again that 'heat' is not a state variable, to which one can 
assign a value in an equilibrium state. In particular, one 
cannot talk of a 'heat content' of the system, despite the 
fact that its existence is strongly suggested by the terms 
'heat current' and 'heat capacity'. For this reason, it is 
preferable to talk about the entropy content of a system, 
its current and its local density, where no such conceptual 
problems exist. 4 



III. DRIFT-DIFFUSION MODEL 

From the point of view of macroscopic thermodynam- 
ics, it is natural to describe non-equilibrium situations 
by spatially varying thermodynamic variables, i.e., the 
local values of the densities and the intensive variables. 
An obvious question is, what is the length scale, below 
which the concepts of a local temperature and a local 
chemical potential are no longer applicable? It is natural 
to identify this length scale with the mean free path A of 
the scattering processes, which are responsible for estab- 
lishing local equilibrium. For smaller distances it is not 
possible to assign a T- or ^-difference. In the simplest 
case, the mean free path is given by 



A 



1 



flscatt C c 



(23) 



where n S catt is the density of scatterers, and a c is the 
cross section of the relevant scattering process. 16 

If the mobile particles in the system move with aver- 
age velocity (\v\) between the scattering events, 1 ' A can 
be translated into a scattering time r via the relation 
A = (\v\) ■ t. The elementary consideration illustrated in 
Fig. 3 shows that the current density associated with a 
balanceable quantity X in linear approximation can be 
written as 



where 



jx = -D-Vx{t,r) , 



D=-{\m 



(24) 
(25) 



Z 




FIG. 3. Elementary derivation of the diffusion constant: 
summing up the four contributions jxz = ±| x(z) ■ (\v\) to 
the z-component jxz of the A"-current density through the 
top and bottom surface of a cube of dimension A, one arrives 
in linear approximation at Eq. 24, if the z-direction is chosen 
parallel to the gradient Vx of the X-density. 



is the diffusion constant. 18 

At face value this derivation appears to be based very 
much on classical physics. Looking more closely, however, 
the only thing really exploited it the condition that the 
quantity X is balanceable, i.e., that it is possible to say 
what is the net amount of X transported in the direction 
of the gradient Vcc(r) of the A-density. For this reason 
the model is extremely robust and holds in both classi- 
cal and quantum physics. Besides on thermodynamics, 
it relies on only two additional concepts, i.e., the trans- 
port velocity (\v\), and the mean free path A between 
scattering events. In quantum physics, 

<MX*) = ^v ge (fc) 

is the group velocity resulting from the dispersion rela- 
tion s(k) of the matter waves, while A is derived from 
the quantum mechanical scattering cross section a c (see 
Eq. 23). 

In order to deal with particle currents, we set X = N, 
and obtain 



j N = -D-Vn(t,r) , 



(26) 



which is known as Fick's 1. law. 

Let us first consider the entropy transport. If we 
choose T and n as independent variables, we can apply 
Eq. 24 to entropy, and obtain the entropy current den- 
sity, and the thermal contribution to the energy current 
density: 



T ■ js 



-T ■ DVs(T, n) 
ds{T,n) 



-T ■ D 



dT 



VT 



ds{T, n) 
dn 



Vn 



-Dc v VT + T 



ds{T,n) 
dn 



3N 



(27) 



where c„ = Tds(T,n)/dT is the heat capacity per vol- 
ume at constant density n = N/V. The first term in 
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Eq. 27 represents Fourier's law, i.e., the conductive ther- 
mal contribution to the energy current with the thermal 
conductivity 



A = Da,, 



1 



nc v (\v\) 2 T 



(28) 



The second term represents the Peltier effect, where at 
constant T the moving particles result in a convective 
thermal contribution to the energy current. The Peltier 
coefficient connecting this contribution and the electric 
current jq = qj^ then reads 19 



n 



T ds(T,ri) 



di 



(29) 



Applying the transport equation to the electric charge, 
we have to choose T and p as independent variables. 
In addition, we have to take into account that even for 
Vn(t, f*) = 0, a so-called drift current can be present, 
which is driven by the electric field. Adding this term to 
the diffusion current (Eq. 26), we obtain for the electric 
current 



j Q = -qD-Vn{t,r)-aV4>(t,>r) 



(30) 



dfi 



dT 



where a is the electric conductivity. In global equilibrium 
T and p, — fi + q<j) are constant. In this case the current 
has to vanish and one obtains the Einstein relation be- 
tween D and a: 



a = -D 



dp, 



(31) 



Note that in equilibrium \i and <j) do not need to van- 
ish separately, a phenomenon, which occurs, e.g., in the 
vicinity of a pn-junction in inhomogeneous semiconduc- 
tors. 

The thermodynamic susceptibility 



dn(T, n) 2 

— 5 = n K T 

op 



(32) 



is closely related to the isothermal compressibility kt = 
— (\/v) ■ dv(T,p)/dp. I propose to call v the particle ca- 
pacity, as it tells us how many particles can be added to 
the system, if the (electro)-chemical potential is raised 
by a certain amount. 20 In systems of charged particles 
particle capacity is related to the electric capacitance. 
In low-dimensional conductors dn(T, [i)jd[i can be very 
small and then provides a contribution to the total in- 
verse capacitance, which cannot be neglected against that 
of the geometrical (electrostatic) capacitance. The sum 
of both contributions determines the ratio between the 
electric charge on a capacitor and the electro-chemical 
potential difference between its electrodes. In the context 
of low-dimensional conductors v is also called 'quantum 
capacitance'. 21 



With the Einstein relation we can rewrite Eq. 30 in the 
more compact form 

h = -f-VM(*,r)-^^§^VT. (33) 

We see that a particle current can not only be driven 
by a gradient of the electro-chemical potential, but also 
by a temperature gradient. The second term in Eq. 33 
describes the Seebeck effect. To identify the prefactor 
in front of VT we employ the Maxwell relation resulting 
from Eqs. 15: 



dn(T,p) _ ds{T,n) _ ds(T,n) dn(T,p) 



dT dfi 
Thus we arrive at 



dn dpi 



(34) 



j Q = -^■Vp(t 7 f)-aSVT(t,f) 

q 



where 



S 



1 ds{T,n) 
q dn 



1 dfi(T,n) 
<? dT 



(35) 



is the Seebeck coefficient, or the thermopower. The sec- 
ond equality in Eq. 35 corresponds to another Maxwell 
equation derived from the free energy density f(T, n) — 
e(T, n)—Ts(T, n). From Eq. 34, and the assumption that 
the same diffusion coefficient D in Eq. 24 holds for both 
the entropy and the particle current density, it follows 
that S and II obey the Kelvin-Onsager relation 22 



n = t ■ s 



(36) 



Applying the general transport equation Eq. 24 to the en- 
ergy density e(T, /i) one can easily verify that the model 
is consistent with Eq. 22, if one uses the homogeneity 
relation Eq. 8 and Eq. 15. 

The main advantage of the drift diffusion model is its 
extreme simplicity and generality. It does not depend on 
the nature of the diffusing particles, and works equally 
well for classical particles, fermions, or bosons. It is also 
independent of the dispersion relation of the particles 
under consideration, e.g., electrons, phonons or photons. 
This means, it can be used for almost all phenomena 
related to the diffusive transport of quasiparticles occur- 
ring in condensed matter physics. Introducing a single 
phenomenological parameter, the diffusion constant D, it 
relates all transport coefficients to thermodynamic sus- 
ceptibilities. 

That single parameter D is also its main deficiency, 
because D usually depends on the energy of the diffus- 
ing particles. As we will see in Sec. VIII, this deficiency 
can be quite easily cured, if particles with different ki- 
netic energy e are associated with different subsystems of 
the gas, having an energy-dependent diffusion constant 
D(e). This additional dependence on energy produces 
in many cases only a prefactor of order unity. Hence, 
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the drift-diffusion model works quite well for the qualita- 
tive consideration of electric, thermal, and thermoelectric 
transport phenomena. 

The drift diffusion model is easily extended to two 
species of particles, e.g., electrons and holes in semi- 
conductors. Another topical example is the transport 
of spin-up and spin-down electrons in the context of 
spintronics, 2 i and spin-caloritronics. 2,25 



IV. CONNECTION TO QUANTUM STATISTICS 

The explanation of the thermodynamic and transport 
properties of matter is the central goal of statistical 
physics. In particular, the thermodynamic derivatives 
entering the transport coefficients in the previous section, 
can be calculated using the methods of statistical ther- 
modynamics. In the following it is shown that the Gibbs 
fundamental form (see Eq. 1), when combined with ele- 
mentary statistical considerations, reproduces in a simple 
way the main results of quantum statistical physics. 

Our starting point is the grand canonical ensemble, 
since this allows an exact implementation of the indistin- 
guishability of quantum particles. This approach allows 
us to use {T, V, fi} as set of independent variables, which 
in the previous section turned out as most appropriate 
for the study of transport phenomena. 

Let us assume that the operators H. and Af represent 
energy and particle number of a given quantum system, 
and have common eigenstates \i) with the eigenvalues Ei 
and iVj. If the state of the system is a statistical mixture, 
in which each eigenstate \i) occurs with the probability 
Wj, the average values of H und Af are given by 26 

E = (u) = Ei Wi ' N = w = E Ni Wi ' 

i i 

(37) 

while the entropy is given by the famous expression 

S = -k B J2w i In Wi. (38) 

i 

Since the probabilities Wi vary continuously, also the av- 
erages E and N form a continuum of real numbers, and 
obey the usual laws of calculus, even if (Af) <C 1. Ther- 
modynamics is now considered as a theory of quantum 
mechanical average values, and the differentials of E, N 
and S read: 

dE = E i d Wi, dN = ^NidWi 

i i 

and dS = -k B J2 ( 1x1 Wi) dWi . 

i 

In the last step we have used the normalization J2i Wi = 
1 of probabilities, which implies Yli = 0- Next, we 
have to determine the probabilities {Wi} such that the 
basis of macroscopic thermodynamics, i.e., the Gibbs fun- 
damental form (Eq. 1) is obeyed. Assuming a constant 



volume V, we write Eq. 1 in the form 



dK({Wi}) = dE -TdS - fidN (39) 
= { Ei + foTlnWi " pNi} dWi = . 



This corresponds to an extremalization of the thermody- 
namic potential if ({Wi}; T, V, p) with respect to the Wj. 
Introducing the normalization condition for the {Wi} via 
the Lagrange multiplier A, we write 



Ei + k B Tln Wi 



fiNi + A} 



dWi = . 



The Wi can now be treated as independent variables, and 
we obtain 



InW 



Ei - fiNi 



A 



knT knT 



(40) 



Hence, we have shown that the probabilities have to fol- 
low the Gibbs distribution 



(41) 



in order to satisfy the Gibbs fundamental form (Eq. 1). 
The quantity 



^(T,y,/2)=exp 



A 



E 



exp 



Ei - fiNj 



(42) 

is the grand canonical partition function. It ensures the 
correct normalization of the probabilities Wi? 7 Z de- 
pends on volume, as the eigenvalues Ei(V) of T-L typically 
vary with V - at least via the the allowed values ki of the 
components of k: ki oc V 2 ^ 3 . 

Next we insert the probabilities into Eq. 38 and obtain 



S = -k B Y, Wi - (- 



Ei - fiNj 
knT 



-lnZ(T,V,ji) 



In Wi 



Using the homogeneity relation Eq. 8 we finally arrive at 
the thermodynamic potential 

K(T,V,fi) = -k B T\nZ{T,V,fi) (43) 
= -p{T,fi)-V . 

Hence the statistical approach allows us to express the 
thermodynamics of a quantum system in terms of its 
grand canonical partition function. The specific prop- 
erties of the system under consideration enter via the 
eigenvalues of its quantum observables % and Af. 



ELEMENTARY FERMI- 
BOSE-SYSTEMS 



AND 



One essential issue of statistical physics is the decom- 
position of many-body systems into simpler thermody- 
namic subsystems, allowing the calculation of K(T, V, p). 
If we restrict ourselves to the simplest case of gases with 
non-interacting quantum particles in a rectangular po- 
tential box of volume V the corresponding Hamiltonian 
reads 



H 



^ e(fc) a\a k , 

k 



(44) 



where the index k represents a wave vector, and Kk is the 
momentum of a particle. The function s(k) is the disper- 
sion relation of these particles and a k ,ak their creation 
and annihilation operators, respectively. The number op- 
erator is simply A/fe = a^ak- When spin is of interest, its 
quantum number a can be simply added to the index fc; 
otherwise it gives rise to an additional factor of 2 in front 
of the sum. 

This Ansatz is very general, because many interacting 
systems can be transformed at least approximately into 
into a Hamiltonian given by Eq. 44. In many cases resid- 
ual interactions between the particles can be taken into 
account as a renormalization of the function e(fc), and a 
finite life time of the quasiparticle states resulting from 
scattering processes. Hence, the Hamiltonian in Eq. 44 is 
relevant for a very broad class of systems, i.e., all systems 
with wave-like excitations. Not only conventional gases, 
but also complex many-body systems such as the lattice 
excitations of solids, Fermi- and Luttinger liquids, su- 
pcrfluids and superconductors, or magnons. Most of the 
quasiparticles of condensed matter physics are at least 
approximately described by the Hamiltonian in Eq. 44. 

The structure of Eq. 44 suggest to decompose an ideal 
gas of particles or quasi-particles into simpler subsys- 
tems. Each subsystem is represented by a single term 
of the sum in Eq. 44, i.e., 

Hk = s(k) a\a k . 

I propose to call these subsystems elementary Fermi- or 
Bose-systems, as they cannot be further decomposed into 
simpler subsystems. In a rectangular potential well 31 the 
elementary Fermi- or Bose-systems share the same vol- 
ume V and have the following properties: 

a) Elementary Fermi-systems have only two eigen- 
states of % and N with the eigenvalues Ei = {0, e}, 
and, respectively, {iVj = 0, 1}. 

The partition function of this system reads accord- 
ing to Eq. 42 



Z F (T,V,p.) = 1 + exp - 



e(fc) — fi 



(45) 
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FIG. 4. Entropy Sk (solid lines) and entropy per particle Sk 
(dashed lines) of elementary Fermi- and Bose-systems with 
characteristic energy e(k). 



b) Elementary Bose-systems have an infinite, but 
countable number of eigenstates of W and M with 
the eigenvalues Ei = {0, e, 2e, 3e, ...}, and Ni = 
{0,1,2,3,...}. 

The partition function is a geometric series in this 
case and reads: 



Z B {T,V,Ji) 



1 



1 — exp 



f e(k)-p, \ 
V k B T J 



(46) 



From the partition functions we can obtain for the aver- 
age particle numbers by differentiating Z k (T, V, ft) with 
respect to /2: 



N k (T,V,fl) 



exp 



V k B T J 



± 1 



(47) 



The N k are the well-known Fermi- (upper sign) and the 
Bose- function (lower sign), respectively. In Fermi sys- 
tems Nk varies continuously between and 1, while in 
Bose-systems varies continuously between and oo. 
For Bose-systems p, has to be always smaller than e(k) 
- otherwise the particle number of the system diverges. 
This divergence of N k at e k = ft is the origin of Bosc- 
Einstcin condensation. It is important to note that the 
N k are average values, around which the particle number 
of the system labelled {k} statistically fluctuates, as op- 
posed to the occupation probability of a single-particle 
state | A:) used in Boltzmann theory. The average par- 
ticle numbers N k are often called distribution functions, 
since they tell how the total number of particles is dis- 
tributed over the different elementary sub-systems. 
The Landau-potential of an elementary Fermi- and 
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Bose systems reads according to Eq. 43 



K k (T,V,jl) = Tk B T ln<M ± exp - 



e(k) — fi 



where the upper sign hold for Fermi-, and the lower sign 
for Bose-systems, respectively. Next we compute the en- 
tropy of the elementary Fermi- and Bose-systems. Using 
Eq. 5 and 43 one finds 



Sk(T,V,fi) 



K k - (s(k)-p)N k 



T 



(49) 



= ±k B ln(l±cxp(-r fc ))± 



Yu 



cxp(r fe ) ± i 



with the abbreviation Y k — (s(k) — £i)/kBT. The upper 
sign refers to Fermi- and the lower sign to Bose-systems, 
respectively. For high energies e ^ fl the second term 
in Eq. 49 dominates, and the entropy per particle ap- 
proaches (see Fig. 4). 



Sfe 



s(k) -/i 



(50) 



The same results are obtained from the thermodynamic 
derivative S k (T, V, fl) = -dK k (T, V, fl) /dT. 

To conclude this section, elementary Fermi- and Bose- 
systems are characterized by two equations of state, 
one for the particle number N e (T, V, fl) (usually termed 
the distribution function), the other for the entropy 
S £ (T,V, fl). The first is ubiquitous in modern physics, 
while the latter is so far not much discussed in present 
literature. This is a pity, as it describes in a very concise 
form the caloric properties of these systems, which are 
the central building blocks of systems with indistinguish- 
able particles. The reason for the lower prominence of S k 
is that the caloric properties of quantum gases are usu- 
ally derived from the energy E e (T, V, p) = £• N e (T, V, p). 
Since e(k) is not a variable, but a characteristic constant 
of the elementary subsystem labelled 'fc', E e appears not 
to be independent from N e , and hence not as an indepen- 
dent equation of state. This is different for the entropy 

S e (T,V,fi). 

In usual terminology the elementary Fermi- and Bose- 
systems are called 'single particle states', as they are 
related to the solutions of the one-particle Schodinger 
equation. The Hilbert space of the many-particle sys- 
tem is represented by the tensor product of the spaces 
of the single particle systems. However, the requirement 
of exchange symmetry resulting from the indistinguisha- 
bility of identical particles eventually prevents the use of 
single-particle systems (whose Hilbert space is spanned 
by the solutions of the one-particle Schrodinger equation) 
as independent building blocks of the many-body sys- 
tem. The exchange symmetry leads to quantum (Fermi- 
or Bose-like) correlations between the single particle sys- 
tems, destroying their statistical independence. These 
correlations are much more conveniently taken into ac- 
count in the framework of second quantization used in 



Eq. 44. The Schrodinger physics still enters in the form 
of the wave functions, which describe the spatial distri- 
bution of elementary Fermi- and Bose-systems. In other 
words, elementary Fermi- and Bose-systems are the ex- 
citation modes of the quantized matter field. These are 
statistically independent, and hence, the thermodynamic 
potential of the gas is just the sum of the thermody- 
namic potentials over all elementary subsystems. The 
particle numbers {N k } and the entropies {S k } of the el- 
ementary subsystems are statistically independent ran- 
dom variables in the sense of thermodynamics, while the 
momenta of different particles are not independent, but 
subject to quantum correlations like the Pauli principle. 

In the next step, we have to assemble the elementary 
Fermi- and Bose-systems to quantum gases. To do this, 
we assume that the fc-vectors compatible with the bound- 
ary conditions lie dense enough in k-space, and convert 
the sum over all elementary systems in into an integral 
over the energies: 29 



? = wrl ddk = Ld J d£9{£) 



(51) 



where g(e) depends on the form of e(k) and the systems's 
dimensionality d; L being the spatial extension of the sys- 
tem in each direction. For spin 1/2 fermions, another fac- 
tor of 2 has to be added or taken into account explicitly, 
if spin phenomena are studied. In usual terminology g{e) 
is called the density of (single particle) states (DoS), from 
our point of view, it is the density of elementary Fermi- 
and Bose-systems on the energy axis. 

The dispersion relation e(k) is (besides the Fermi- or 
Bose character) the main characteristic of the combined 
system, i.e., the specific quasi-particle gas under investi- 
gation. It is the only feature of our description that has a 
counterpart in classical physics (i.e., the function E(P)). 
The difference between the classical and the quantum 
point of view is the following: rather than saying that the 
same (classically distinguishable) particle is accelerated 
by external forces or a scattering process with another 
particle, the quantum point of view is that particles in 
one elementary Fermi- or Bose-System are annihilated, 
while a (possibly different) number of particles is gener- 
ated in other elementary Fermi-or Bose-System (of course 
under the constraint of the applicable conservation laws 
for the transitions). The dynamical variables are not any- 
more the positions and momenta of individual particles 
(which do not exist because of indistinguishability) , but 
the particle number N kl the entropy S k , and the energy 
E k = e(k)N k of the elementary Fermi- and Bose-systems. 

Using the integral representation of Eq. 51, it is 
straightforward to compute the thermodynamic poten- 
tial of the quantum gases from that of the elementary 
subsystems (see Eq. 48): 

K(T,V,p) =-k B TV deg(e)K £ (T,fi) , 
Jo 

where V — L d . From the experimental point of view, it 
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is more useful to compute the thermal 

/>oo 

n(T,fi) = / dsg(e)-N E (T,p,) 
Jo 



(52) 



and the caloric 
e(T,/2) 



deg(e)- N e (T,p,)-e (53) 



equation of state. As a variant of the caloric EoS, the 
entropy density of an ideal gas is at all temperatures 
given by 



s(T,p) 



deg(e)-S £ (T,ft) 



(54) 



Interestingly, all the entropy of the gas comes from the 
statistical mixture of states with different particle num- 
bers within the same mode (or the same elementary sub- 
system). The combination of the elementary systems 
corresponding to different /c-directions to the composite 
system does not enhance the entropy of the gas beyond 
the sum of the entropies of its elementary subsystems. 
This means that the entropy of a gas results exclusively 
from the statistical mixture of states with different parti- 
cle numbers within one elementary subsystem! This will 
be important for the following section, where we will see 
that the reduced dimensionality of nanostructures leaves 
the properties of the individual elementary subsystem un- 
changed. 

From the above equations of state all the thermody- 
namic susceptibilities entering the transport coefficients 
of the drift-diffusion model (see Sec. Ill) can be com- 
puted in equilibrium, where all these systems share the 
same temperature and chemical potential. For degener- 
ate electrons, e.g. in three dimensions, one obtains in 
effective mass approximation the DoS 



with the Fermi energy 



3 n 



,1/3 



2e F {n) 



ep(n) = ji{T = 0, n) = 



ft 2 (37r 2 n) 2 / 3 
2m 



and one arrives at the well known thermodynamic sus- 
ceptibilities: 



2 dn{T,ii) 

n k t = = g(s F ) 

o\x 

de{T,n) ds(T,n) 
c v {l,n) = — — — = 1 - 



(55) 



dT 
7T 2 nk 2 B T 



dT 



2 £p(n) 
ds(T,n) _ ir 2 nk 2 B T _ s(T,n) 
dn 6 £p(n) 3 



s(T,n) (56) 
(57) 



where kt the isothermal compressibility, and c v the ther- 
mal capacitance per volume at constant particle density. 



On the other hand, for dilute gases, e.g., the electrons 
in semiconductors one obtains 



n 2 kt — 


dn{T,n) 
dfx 


c v 


de(T, n) 
dT 


s(T,n) = 


nks |ln ( 


ds(T, n) 
dn 


s— k B = 



k R T 



T- 



dT 
]T 3 ' 2 \ 6 



2/ ' 
(ill 



(58) 
nk B (59) 
(60) 



,(61) 



where j = 2 • (™fcs/27r/i 2 ) 3 / 2 is called the chemical con- 
stant of the electron gas (the prefactor 2 comes from 
spin). Eq. 60 is also known as the Sackur-Tetrode equa- 
tion. Eq. 61 provides a reasonable estimate for the ther- 
mopower of a doped (non-degenerate) semiconductor in 
the T-regime, where most dopants are ionized. 32 Because 
of the small factor ksTp/sF in Eq. 57 the thermopower 
of metals is much smaller than that of semiconductors. 

The fact that the state of macroscopic system with 
~ 10 23 internal degrees of freedom can be specified by 
only three independent variables (here T, V, and fj,) is a 
consequence of the thermal and chemical equilibrium be- 
tween the elementary subsystems. The composite char- 
acter of macroscopic systems becomes visible, if the equi- 
librium between the subsystems is disturbed, e.g., by a 
laser, which selectively increases the population of ele- 
mentary subsystems at higher energy at the expense of 
those at lower energies. The equilibrium is restored by 
inelastic recombination processes. 

The question is whether the elementary Fermi- and 
Bose systems are fundamentally relevant as separate en- 
tities, or are just a technicality, which allows an easy 
computation of the EoS. This question is equivalent to 
the question of why the grand canonical approach to sta- 
tistical thermodynamics should be preferred to the micro- 
canonical and canonical ones. In equilibrium and in the 
thermodynamic limit, where the fc-space can be consid- 
ered as a three-dimensional continuum, also the micro- 
canonical and the canonical approaches to thermodynam- 
ics work reasonably well. This changes when going be- 
yond these restrictions as discussed in the next two sec- 
tions. Historically, the characteristic statistical fluctua- 
tions of energy and particles numbers in these ensembles 
strongly biased the common view. Before the advent of 
quantum mechanics fluctuations were considered possi- 
ble only for open systems. Statistical fluctuations were 
a nuisance within classical physics, and their appearance 
within the canonical and grand canonical ensemble was 
viewed as an artifact, which ought to be removed by 
the thermodynamic limit, i.e., the limit V, N — > oo at 
constant particle density n. From this perspective, the 
micro-canonical approach is often considered as the most 
fundamental one, leading to the perception that thermo- 
dynamics as a whole works only within this limit. In 
quantum physics, however, statistical fluctuations con- 
stitute an unavoidable element of physics (see section IX 
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for further discussion). 



VI. BALLISTIC QUANTUM TRANSPORT OF 
PARTICLES 



In this section, it is shown that the elementary Fermi- 
and Bose-systems (which may look artificial at first sight) 
are very useful to understand transport properties in 
reduced dimensions, where the drift-diffusion model is 
entirely inapplicable. This is the regime of mesoscopic 
transport, 33 where the size of the conducting object is 
smaller than the mean free path A. More precisely, it 
is even sufficient, if the inelastic mean free path Aj„(T), 
corresponding to those scattering processes resulting in 
the dissipation of energy, is larger than the sample size. 

Elastic scattering, i.e., scattering at constant energy, 
cannot create entropy because any increase AS of en- 
tropy requires the amount AE = T AS of energy, which 
by definition is not available. Instead, the elastic scat- 
tering modifies the underlying wave function, i.e., it 
changes the spatial distribution of the elementary Fermi- 
or Bose systems. The plane wave functions of ideal 
gases, or Bloch-wave functions of crystals, respectively, 
are then replaced by complex interference patters, which 
are known as speckle patterns from laser physics. 

At low temperatures T < 1 K the inelastic mean free 
path Ai n (T) is typically in the micron regime, while in 
high mobility semiconductor also the elastic mean free 
path easily exceeds the micron range, and enables the 
study of truly ballistic transport. In this regime it is 
possible to experimentally realize some simple textbook 
quantum systems, where plane waves are scattered off 
tunable potential barriers. 

The simplest case to consider is that of a quantum 
wire: a quantum wire can be realized for photons (wave 
guides or optical fibers comparable with the wave length 
of electromagnetic waves), for phonons (narrow beams 
with a diameter comparable to the wavelength of ther- 
mally excited phonons), and for electrons (semiconduc- 
tor heterostructures or carbon nanotubes with a diam- 
eter comparable to the Fermi wavelength), to list a few 
examples. The common element in these examples is 
the fact that the set of allowed fc-vectors forms not a 
three-dimensional continuum anymore, but consists of 
one, or a few one-dimensional sub-continua, which prop- 
agate particles, entropy, and energy along the wire (say, 
in x-direction), while the transverse part of the wave 
function is discrete, resulting from the strong confine- 
ment of the system in the transverse direction. The one- 
dimensional sub-continua of elementary Fermi- or Bose- 
systems are also called transport channels. For simplicity, 
let us first assume that we have only one transport chan- 
nel, as sketched in Fig. 7. Initially, we also assume that 
the quantum wire is perfectly transmitting (T(e) = 1): 
smaller transmission coefficients T(e) < 1 are easy to 
take into account in the next step. 

The transport of photons, phonons or electrons can 



then be viewed as a scattering problem: a beam of parti- 
cles, emanating from two particle reservoirs connected to 
the left and right end of the wire, is either transmitted or 
reflected back. The elementary Fermi- or Bose-systems 
in the wire break up in two subsystems: right movers 
and left movers, which propagate particles, entropy, and 
energy with the dynamical velocity 



v(k) 



de{k) 
hdk 



(62) 



and are populated according to the temperature and the 
electrochemical potential of the left and right reservoir: 34 

Ne;L,R = N e (T L<R , p, LjR ) . 

In this approach the elementary Fermi- and Bose-systems 
can be visualized as conveyor belts for energy, entropy, 
particles, momentum, and spin, which transport these 
quantities ballistically, until an elastic, or inelastic scat- 
tering event occurs. 

If the two reservoirs differ in temperature or in 
(electro)chemical potential net currents of E, S, and N 
will flow. Despite being in a non-equilibrium state as a 
whole, the currents flowing through the wire of length L 
are perfectly described by the thermodynamic properties 
of the two (left- and right-moving) subsystems. Following 
Landauer and Biittikcr 30 ' 36 we can write 



In = 



1 

L 



(63) 



.fc>o 



k<0 



de g(e) v(e) {N £ (T L , fi L ) - N £ (T R , ji a )) 



where v(k) = — tf( — kj. In the second step we have 
evaluated the sum in a continuum approximation using 
the one-dimensional DoS g(e) = (1/tt) ■ dk x (e)/de for 
propagating modes in one dimension. 

If we plug this into Eq. 63 the energy dependence of 
the DoS and the velocity cancel and we obtain the sur- 
prisingly universal result 



/ \ - ^ I ds(N s (T L ,ji L )-N £ (T R ,ji R )) , (64) 



which is valid for both Fermi- and Bose-systems and in- 
dependent of the functional form of the dispersion re- 
lation e(k). In solid state nanostructures, phonon and 
photon currents are usually not detected by measuring 
electric current or counting particles, but as a thermal 
(i.e. 'heat') current. We leave the treatment of ther- 
mal currents to the next section and specialize now to 
charged systems, in order to compute the electric con- 
ductance and the thermopower of quantum wires. 

If the wire hosts several transport channels with e- 
dependent transparencies T n (e), the charge current Iq — 
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qljy assumes the more general form 



J2_ 
kH 



J deT n (e) (N s (T L ,fl L ) - N e (T R ,fl R )) 



(65) 

Here, we limit ourselves now to the linear response 
regime, and assume that the wire is symmetrically bi- 
ased, i.e., 



N e -,L,R= N £ (T ± AT/2, n ± qU/2) , 



(66) 



where AT and U are the applied temperature and volt- 
age bias, and the upper and lower sign refer to the left 
and right reservoir, respectively. Then we can Taylor- 
expand the N E . l r of the reservoirs around the averages 
(Tj, + Tr)/2 and (fix, + jl R )/2 and obtain in linear ap- 
proximation (fcflAT, qU <C ksT): 



N s (T L ,fi L ) - N e (T R ,fi R ) 
dN(Y) 1 



dY k B T 



dN(Y) 
dY 

qU 



AY 



(67) 



T 



-AT 



AY 



where Y — (e — /i)/fcflT, and 

dN(Y) exp(r) 



dY (exp(y)±l) 
The charge current then reads 

dN(e) t TT e-fj, 

o 



(68) 



OO 

J Q = ^E J d£ T n (e) 



de 



qU - '-^AT 



(69) 

In case of fermions, a Sommerfeld expansion of the inte- 
gral in Eq. 69 leads to 



where 



In = G-U + GS-AT 



G = G -2j2r n (^) , 



(the factor 2 takes into account spin degeneracy) and 



G n = 



38.74 nS 



1 



25.8 kfi 



(70) 



is the universal conductance quantum. This is a seminal 
result of mesoscopic physics found experimentally first in 
gate-defined quantum point contacts. ' 38 

The thermoelectric counterpart of the conductance 
quantization, i.e., the Seebeck coefficient of a quantum 
wire, or quantum point contact is given by 



S = 



^ • 2 > 

3 qh ^ 



d(]nT n (s)) 



de 



(71) 



Also this result has been experimentally confirmed first 
in quantum point contacts. 39 



VII. BALLISTIC QUANTUM TRANSPORT OF 
ENTROPY 

The first efforts to transfer the ideas of ballistic elec- 
tron transport to thermal transport originate from Imry 
and Sivan 40 and Butcher. 41 To describe the thermal 
transport through quantum wires, we can write down 
an expression for the entropy current that is the ther- 
mal analogue of Eq. 64, but contains the entropy Sk- l,r 
(Eq. 49) of the elementary Fermi- or Bose-systems rather 
than their particle numbers N e -L, R - Assuming that the 
entropy propagates in each elementary Fermi-, or Bose- 
system at the same velocity v(k) as the particles, the 
entropy current reads: 



Is = 



OO 

^ J de (S e (T L ,p, L ) - S e (T R ,p R j) 



(72) 



For the same symmetric bias (see Eq. 66) one finds in 
linear approximation 



S E (T L ,fl L )-S s (T R ,fi R ) ^^pAY 

= ^^(qU-^AT 
dY k B T V T 

where again Y = (e — ii)/k R T, and 

dS(Y) t v exp(y) 



(73) 



dY 



k R Y 



(exp(y)±l) 



(74) 



The upper sign holds for fermions, and the lower for 
bosons. Interestingly, this result differs from Eq. 68 only 
by the extra factor ksY. For the entropy current we then 
find in linear response 



Is 



de T„(e) 



H dN(e) 



T 



de 



W-^at 



(75) 

The first term (oc Afi) in this equation is driven by the 
voltage bias, and constitutes the ballistic analogue to the 
Peltier current in Eq. 27. The second term (oc AT) de- 
scribes the entropy current driven by the T-difference. 

The same result is obtained, if one extends the deriva- 
tion of In to Ie and computes Is via Eq. 22: 

Is = ^(Ie-iUn) = ^{n- U + C-At) (76) 

Note that the identification of Ig with (Ie—JIIn)/T holds 
only in the linear response regime. 

For fermions, we can again evaluate the integrals in 
Eq. 75 within the Sommerfeld approximation, and obtain 
for the Peltier coefficient 



9 „ 



de 



e—^i 
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and for the thermal conductance 

C = T-2L ^T„(m) , 

n 

where the prefactors 2 accounts again for spin- 
degeneracy. The quantity Lq is the entropy conductance 
quantum corresponding to 

L = ~ f = 0.9456 pW/K 2 , (77) 

implying that in ballistic quantum wires the entropy con- 
ductance L — L/T is quantized in units of Lq. Compared 
to the quantum of electric conductance, q 2 is replaced by 
(nk B f/3. 

As discussed in the context of the drift-diffusion model 
in Sec. Ill the Seebeck- and Peltier-coefficients are con- 
nected by the Kelvin-Onsager relation 

n = T-S. (78) 

This results from the Maxwell-relation 
dN £ (T,fl) = dS £ (T,p) = exp(V) 

dT dp T(exp(F)±l) 2 ' 1 ' 

between the derivatives of the two EoS (Eqs. 52 and 54), 
and the fact that the transmission coefficients T(e) de- 
termine all transport quantities. 

These results are again very general, as they de- 
pend only on the energy dependence of the transmis- 
sion coefficients, but neither on the dispersion relation 
nor on the particle statistics. After some more quali- 
tative experiments 44 a quantitative experimental inves- 
tigation of the thermal conductance of quantum point 
contacts has been performed only recently, exploiting 
the thermopower of quantum point contacts for lo- 
cal thermometry. 45 Very recently, also chiral thermal 
transport in the integer quantum Hall regime has been 
demonstrated. 46 

The so far most relevant cases for bosons deal with 
phonons and photons. In this case we can set /i = 0. 
The evaluation of the corresponding Bose-integral re- 
sults in a quantized thermal conductance with the very 
same entropy conductance quantum L as in the case of 
fermions. ' The case of the quantized entropy conduc- 
tance by phonons was first addressed by the beautiful ex- 
periments by Schwab et al. 49 A few years later, Meschke 
et al. considered the case of microwave photons. 50 The 
latter case is of particular importance for instrumentation 
in mesoscopic physics. First, it explains why a careful 
filtering of the measurement leads at low temperatures 
is required: in a cryogenic setup the wires bring down 
energy and entropy not only via electronic and phononic 
thermal conduction, but also via thermal photons. These 
photons may not carry enough energy to heat macro- 
scopic objects like the thermometer, but they can still 
induce jumps of charge carriers in electronic traps, or re- 
lease single electrons from a quantum dot. Second, the 
photon case is interesting, because the techniques of mi- 
crowave engineering provide possibilities to manipulate 
the photon thermal transport at the mesoscopic scale. 51 
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FIG. 5. The left- (right-) propagating elementary 

Fermi- or Bose-systems labelled 'fc' emanating from the 
right (left) volume element of size contribute an amount 
VkXk{TL,R, P-l,r) to the total X-current density jx- 

VIII. CONNECTION TO DIFFUSIVE 
TRANSPORT 

The concept of elementary Fermi- and Bose-systems 
can also be used to analyze the case of diffusive trans- 
port more accurately than in the drift-diffusion model. 
To do this, we consider in analogy with section III two 
adjacent volume elements of linear dimension A = \v\ ■ r, 
which arc illustrated in Fig. 5, and assume that an inelas- 
tic scattering mechanism exists, which establishes local 
thermal and electro-chemical equilibrium on the scale of 
A. Volume elements of this size represent the smallest 
possible units, which can be assumed to be in local equi- 
librium. For a given T- or /i-gradient, two adjacent cubes 
exhibit the minimal T- or /2-difference possible. 

Between collisions the elementary Fermi- and Bose- 
systems propagate energy, entropy and particles ballisti- 
cally, as they do in the one-dimensional ballistic quantum 
wires discussed in the preceding sections. To account for 
the three dimensions we have to average over the different 
fc-directions, in order to determine the current density of 
any balanceable quantity A in a given direction. 

Similar to section III, we write the A-current density 
through the interface between the two cubes as the dif- 
ference 

jx = ^2v(k) ■ (x k (T L ,p L ) - x k (T R ,p R )\ (80) 

k 

of left and right propagating currents, where x k is the 
contribution of each elementary subsystem to the total x- 
density depending on the values of T^ R and Pl,r within 
two adjacent cubes of size (see Fig. 5). The simi- 
larity of this expression to those used in ballistic trans- 
port is not accidental, but results from the fact that 
the propagation of any balanceable quantity A is bal- 
listic over distances smaller than the mean free path A k . 
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In the examples considered here Xk is the contribution 
rife = Nk(T,p,)/V or Sfe = Sk(T,p,)/V of an elementary 
Fermi- or Bose-system with the wave vector k to the par- 
ticle density, or the entropy density, respectively. 

In contrast to the more elementary treatment of the 
diffusive limit in section III we now take into account 
that Afc = v(k) ■ Tj, depends via both Tfe (see Ref. 52) and 
the velocity v(k) of propagation (see Eq. 62) between 
collisions depend on e, or more precisely on k. In linear 
approximation, the difference Ax k reads in analogy to 
the preceding sections VI and VII 



%k,L — x k,B, 



dx k (Y) 
dY 



AY 



(81) 



where again Y = (e — fij/keT, and AY is the variation 
of Y over the distance of ballistic propagation between 
scattering events. Similar to Eq. 67 AY is given by 



ay, = — (vn- 5 ^ 

k B T 1 T 



VT) -v(k)T(k) . (82) 

A fc 



Here, |Afe| is the mean free path associated with the scat- 
tering of quasiparticles in the elementary subsystem with 
wave vector k. 

Plugging AYfe into Eq. 81 we obtain in the continuum 
limit 



jx = / de g{e)D(e) 



dx(e) 
de 



where the diffusivity tensor 



V/2- 



S — fJL 

T 



VT , (83) 



D(e) = JJ d9dipv(e,6,ip)®v{e,6,<p)-T(s,6,(p) (84) 

is the energy-dependent tensorial equivalent of the dif- 
fusion constant (averaged over the angles 9 and tp on 
a sheet of constant e(k)), which takes into account the 
shape and angle anisotropy of the dispersion relation e(k) 
and the scattering time r(k). 30 For an isotropic e(k) and 
r(fc) Eq. 84 reduces to Eq. 25. 

Equation 83 constitutes the linearized transport equa- 
tion for any balanceable quantity X, and is in per- 
fect agreement with the result from the solutions of the 
Boltzmann equation in relaxation-time approximation in 
textbooks. 30 However, it is remarkable that its derivation 
here is based not on the classical concepts of trajectories, 
but on the sole assumption of the existence of a scatter- 
ing mechanism ensuring that the elementary Fermi- and 
Bose systems propagate (on average) ballistically from 
opposite faces through a cube of size A 3 have values of 
Xk(T,p) according to the local values of T(f), fi(r), and 
T(r + A fe ), p,(f + Afe), respectively. 

The evaluation of Eq. 83 for the particle and the en- 
tropy current reproduces the Drude formulas for the elec- 
tric (Eq. 31) and the thermal (Eq. 28) conductivity, as 
well as the Wiedemann-Franz law. The evaluation of 



the thermopower and the Peltier coefficient within the 
Sommerfeld expansion leads to a slightly modified result, 
when comparing to the drift-diffusion model (see Eq. 35). 
While the thermodynamic derivative dn{T,[i)/dT enter- 
ing Eq. 35 involves only the derivative dg(e)/de (because 
the diffusion constant D is assumed to be independent 
of £■), the evaluation of Eq. 83 contains the derivative 
d{g{e)D(e)) / de . The thermopower is then given by the 
Mott-formula 



k%T d(ln[g(e)D(e)}) 



de 



(85) 



which is the diffusive analog of Eq. 71. If both g(e) oc 
e a and D(e) oc obey a power law, the logarithmic 
derivative in Eq. 85 assumes the value 



d( In g (e)D(e)) 



de 



■0 



E=fl 



as opposed to the result a/ep obtained within the drift- 
diffusion model. For free electrons in three dimensions we 
have a = 1/2 and j3 = 1, if we assume the scattering time 
r to be independent of e. In this simplest approximation, 
the result of Eq. 85 is a factor of three larger than Eq. 57, 
and wc have the curious relation 



s(T, n) 



(86) 



This result is sometimes exploited to estimate the mo- 
lar entropy of strongly correlated Fermi systems from 
measurements of the thermopower. In doing so, one has 
of course to be aware, that Eq. 86 heavily relies on a 
quadratic e(k) and the energy-independence of r. It is 
certainly much safer to extract s from measurements of 
the molar heat capacity. 

So far we have restricted our consideration to an iso- 
lated degenerate Fermi gas. A more accurate modelling 
has to take into account that a solid usually contains 
several different Fermi- or Bose-gases. Besides the elec- 
trons, there are phonons, and (in magnetic solids) also 
magnons, or localized magnetic moments. The interac- 
tion between the different quasiparticle systems provides 
usually scattering mechanisms with an energy-dependent 
scattering time r(e), which may deviate from a power 
law, or at least, affect the value of /3. In addition, 
there can be drag phenomena, like electron-phonon, and 
electron-magnon drag, which can substantially compli- 
cate the behavior of S. 53 

The mean free path Afe (given by Eq. 23) can be calcu- 
lated from the solution of the quantum mechanical scat- 
tering problem of the relevant particles. It is quite inge- 
nious that the model describes the irreversible process of 
conduction, without the need to explain how entropy is 
actually generated in the scattering process. As we will 
discuss in section IX D, the scattering cross section is cal- 
culated within standard Hamiltonian quantum theory, in 



15 



which the time evolution is always reversible. Neverthe- 
less, as the scattering enters only via a scattering proba- 
bility (and not via a probability amplitude) in Eq. 83, it is 
assumed that the phase coherence in the scattering pro- 
cess is quenched. This looks like a rather arbitrary trun- 
cation of the Hamiltonian time evolution that has some 
resemblance with the quantum mechanical measurement 
process. It is this (here manually imposed) quenching 
of the phase coherence, which gives rise to the 'classical' 
character of the Boltzmann-like transport theory. In the 
next section it is discussed, how phase coherence affects 
the transport. It turns out experimentally, that it is the 
inelastic scattering processes, which account phenomeno- 
logically for the loss of phase coherence. 



IX. DISCUSSION 

A. Applicability of Thermodynamics 

It is a widely spread opinion that the concepts of ther- 
modynamics are applicable only in equilibrium, and for 
large systems with many microscopic degrees of freedom. 
The limit of large systems is also called the thermody- 
namic limit, where N, V — > oo, and n — N/V — const. 
These restrictions result from the custom to define en- 
tropy by Boltzmann's famous formula 

S = k B \nQ(E,V,N) , (87) 

where O is the number of microstates accessible to N 
particles in a volume V with a given value of the to- 
tal energy E. This strategy (when properly combined 
with the principle of indistinguishability) indeed allows 
to calculate the function S(E, V, N), which is equivalent 
to E(S, V, N), and hence constitutes a starting point for 
equilibrium thermodynamics in entropy representation. 
However, in so-called open systems, which exchange en- 
ergy and particles with the environment, the number f2 
of microstates cannot be defined anymore, and hence it 
seems that thermodynamics is inapplicable to transport 
situations. 

On the other hand, the presentation of the preceding 
sections shows that thermodynamics can well be applied 
to a large variety of non-equilibrium and transport situ- 
ations, provided that a proper decomposition into simpler 
sub-systems is chosen. Each of these elementary sub- 
systems is at least approximately in thermodynamic equi- 
librium with one of the reservoirs, while the thermal and 
chemical equilibrium with its fellow sub-systems propa- 
gating in the opposite direction can be disturbed. 

To be more specific, in a ballistic quantum wire such 
as the one-dimensional subbands formed in high mobil- 
ity two-dimensional electron systems (see Fig. 6) , the set 
of right- (lcft-)moving elementary Fermi, or Bose sys- 
tems are in equilibrium with each other and the left 
(right) reservoir, but the equilibrium between left and 
right movers is disturbed by the applied bias voltage, or 



the temperature difference, respectively. The same holds 
locally in a macroscopic piece of matter subjected to a T- 
or /2-gradient: also here left- and right moving elementary 
Fermi- and Bose-systems are out of equilibrium, but the 
amount of entropy and particles stored in the elementary 
Fermi- or Bose-systems propagating in the same direction 
from a certain point in space is to a very good approxi- 
mation given by the local values of the temperature and 
(electro)-chemical potential. 

This explains, why the current densities of entropy and 
particles are determined by the diffusion coefficient and 
the local (equilibrium) values of the derivatives of densi- 
ties nk(T,p,) and Sfc(T,/z), respectively. 

In this sense the thermodynamic concepts retain their 
relevance, irrespective of a reduced dimensionality or the 
absence of global equilibrium. This is also intuitively 
clear, as we don't question the validity of the concept of 
water temperature based on the undisputable absence of 
thermal equilibrium between the Mediterranean and the 
polar sea. That the same concepts remain valid in quan- 
tum wires with perfect transmission down to the atomic 
scale, e.g., in highly transparent atomic point contacts, 
appears more surprising, but has been demonstrated in 
a variety of beautiful experiments. 57 

It has to be pointed out that a large body of literature 
considers the Hamiltonian dynamics of small quantum 
systems coupled to reservoirs, striving for a modelling 
of nanoscale thermodynamic processes (see, e.g., Ref. 58 
and the references therein) , including the transfer of work 
and heat between nanosystems and the reservoirs. In 
contrast the present work is focused on simple transport 
phenomena, and tries to elucidate the conceptual foun- 
dation, on which the standard approach to both ballistic 
and diffusive transport is so successful. 



B. Breakdown of Local Equilibrium 

The Landauer-Buttiker approach to transport restricts 
itself to the particular case of a nanoscale constriction 
between macroscopic reservoirs, in which the thermody- 
namics of the elementary Fermi- or Bose-systems in the 
constriction is governed by the reservoirs, and can be 
separated from the transmission properties of the con- 
striction expressed by the set {7^(e)} of transmission 
coefficients. 59 This separation holds best in the linear 
regime, and as long as energy and particle number of 
the elementary subsystems differing in the characteristic 
energy e are statistically independent, because no inelas- 
tic scattering induces an exchange of energy, entropy and 
particles between them. In this way the population of the 
reservoirs, which are by definition in thermal equilibrium, 
is transferred to the elementary subsystems. These can 
be considered in thermodynamic equilibrium with one of 
the reservoirs, while there is no equilibrium between ele- 
mentary subsystems charged by different reservoirs. 

Under these conditions, the propagation of energy, en- 
tropy and quasiparticles within the constriction is gov- 




FIG. 6. Schematic of single channel quantum wire connected 
to two reservoirs for energy, entropy, and particles. The one- 
dimensional character of the transport can be realized by the 
formation of edge states in a quantizing magnetic field, i.e., in 
the quantum Hall regime. In this case the magnetic provides 
also a spatial separation between left- and right-movers. The 
elementary Fermi-or Bose systems are charged with energy, 
entropy, and particles via the left (red) and the right (green) 
reservoir, respectively. They are in thermodynamic equilib- 
rium with their source reservoirs, but not with each other. 

erned by reversible Hamiltonian dynamics, while the con- 
ceptionally difficult irreversible equilibration of the in- 
jected quasiparticles within the reservoirs is irrelevant 
for the transport properties of the constriction! This ex- 
plains the success of the Hamiltonian dynamics in the 
description of this category of transport processes. It is 
interesting that this success can (within the semi-classical 
approximation) also be transferred to the macroscopic 
case, where reservoirs are absent, and a sufficiently strong 
inelastic scattering ensures local thermodynamic equilib- 
rium, if one averages over volumes larger than the mean 
free path (see Sec. VIII). 

If elastic scattering occurs within the phase coherent 
region, the elementary subsystems labelled {k} in the 
absence of scattering have to be replaced by more com- 
plicated ones, which in usual terminology are called scat- 
tering states, and consist of the incoming and the two (or 
more) outgoing waves. Such a complex wave pattern still 
forms one elementary Fermi- or Bose system that cannot 
be decomposed further. In particular, it is impossible to 
provide a local thermodynamic description of the quan- 
tum wires left and right of the QPC in Fig. 7. The QPC 
partitions the flux of energy, entropy and quasi-particlcs 
emanating from each reservoir according to the transmis- 
sion coefficient T{e) into the two corresponding outgoing 
fluxes. As the particles emanating from different reser- 
voirs are incoherent, and thus statistically independent, 
their average particle numbers in the elementary Fermi- 
and Bose systems flowing into the right reservoir simply 
add up according to 

ivp = T(e)-N^(T L ,Jl L ) + (l-T(s))-Nll(T R ,fi R ) , 

(88) 

while we have for the left reservoir 

N'lT = (l-T(e))-N^(T L ,fi L ) + T(e)-N^ k (T R ,fi R ) . 

(89) 

Here N^ q denotes the equilibrium particle numbers in 
the reservoirs. According to Eq. 64 the non-equilibrium 
particle numbers N^ cq determine the current. Analogous 
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FIG. 7. The same quantum wire as in Fig. 6, but interrupted 
by a quantum point contact (QPC) with transmission coef- 
ficient T(e). The scattering states of the particles emanat- 
ing from the left (red) and the right (green) reservoir, re- 
spectively, have to be considered as one elementary Fermi- 
or Bose-systems, which cannot be further decomposed into 
subsystems. 



expressions hold for Sk, and all other balanceable 
quantities of the system. 

Being a superposition of equilibrium particle numbers, 
and entropies with different {T, /!}, the particle numbers 
Nk in the outgoing branches do not satisfy Eq. 47 any- 
more. This signals the break-down of local equilibrium. 
Similarly, the entropies T{e)Sk and (l — T(e))Sk pro- 
pagated by the outgoing partial waves after the scatterer 
are lower than the entropy of an equilibrium mode with 
the same temperature of the left reservoir (see Eq. 49) 
and particle numbers T{e)Nk and (l — T(e))Nk, respec- 
tively. 

One may be tempted to consider a single quantum wire 
between the scattering region and one reservoir as a 'sys- 
tem' on its own, and ask for its thermodynamic quan- 
tities, e.g., its particle number or entropy. This object, 
however, is certainly not an independent sub-system in 
the sense of thermodynamics, because its particle num- 
bers and other physical quantities are correlated with 
those of the other quantum wires connected to the con- 
striction. The correlation of the particle numbers can be 
experimentally accessed, e.g., as an anti-correlation be- 
tween the particle currents in the two reservoirs. 61 Thus 
it has to be kept in mind that the transport modes, i.e., 
the set of elementary Fermi- and Bose-systems hosted by 
the scattering region and the quantum wires cannot be 
decomposed in smaller subsystems as long as their phase 
coherence is not destroyed by inelastic processes. 

A particularly nice experimental demonstration of such 
out-of-equilibrium physics was recently obtained using 
the one-dimensional edge channels in the quantum Hall 
regime. 62,63 There it was possible to measure the e- 
dependence of the particle numbers N £ at different dis- 
tances (0.8-30 /xm) from a quantum point contact. Close 
to the quantum point contact a double step shape of N e 
was observed, which resulted from the superposition of 
the transmitted and reflected electrons emanating from 
of reservoirs at different electrochemical potential (see 
Eqs. 88 and 89). At larger distances inelastic scattering 
processes restored the local thermal and chemical equi- 
librium, and N e behaved Fermi-like. The measured lo- 
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cal temperature was higher than that in the reservoirs, 
but lower than expected from energy conservation argu- 
ments. One possibility to explain the apparent loss of 
energy are additional neutral modes, which are predicted 
to exist, when two (spin-degenerated) edge channels are 
present. 64 

A breakdown of local equilibrium can also happen in 
the diffusive regime. As it was shown in a series of ex- 
periments on diffusive wires of length L in the /im-regime 
this is possible at low temperatures, if the average diffu- 
sion time td = L 2 /D through the wire is shorter than 
the time required for energy relaxation in the wire. 65 This 
means the energy relaxation occurs predominantly by in- 
elastic scattering in the two macroscopic reservoirs serv- 
ing as contacts for the wire. If the energy relaxation in 
the wire is entirely negligible, no particle exchange occurs 
between elementary Fermi-systems of different character- 
istic energy e - i.e., elastic scattering processes dominate 
in the wire. The effect of the elastic scattering is an effi- 
cient randomization 66 of the momentum distribution in 
the wire, which corresponds to intense exchange of par- 
ticles (and entropy) between elementary Fermi-systems 
with the same e, but differing in fc-dircction. 6 ' In absence 
of inelastic scattering the elementary Fermi-systems with 
the same e are not in local equilibrium, because N e and 
S £ at any point in the wire have contributions from both 
reservoirs. 

To quantify these considerations, it is simplest to con- 
sider the diffusion equation for the particle density n £ {r) 
in elementary Fermi-systems with the same e (r is the po- 
sition along the wire) with a source term S x e describing 
the exchange of particles with elementary Fermi-systems 
with other e. The diffusion equation results from the 
combination of the continuity equation Eq. 20 with Fick's 
law (Eq. 26) for the particle densities n £ (r) in the com- 
posite sub-system containing all elementary Fermi sys- 
tems with the same characteristic energy e: 

dn e {rt) +di ^ = s _ (9Q) 

ot 

The source term Ejv> is called the scattering integral, and 
consists of a sum over all possible transitions between 
one elementary Fermi-system and the others with the 
quantum-mechanical transition probabilities W £ . £ '. 65 For 
isotropic scattering only the energy dependence of the 
W £t£ > is relevant, and it reads 

Ejv e = - J de' g{e)g(e') x 

{w s , e m e (l - r V ) - W^ e rv(l - n £ )} . (91) 

In this case a relaxation time approximation is not suf- 
ficient to describe the experimental data, which are ac- 
curate enough to extract the energy dependence of the 
scattering probabilities W £ . e '. 

If inelastic scattering occurs only in the reservoirs, but 
is negligible in the wire, Sjv £ can be neglected and the 



resulting n £ is a two-step function composed of the two 
Fermi functions of the reservoirs with a weight factor that 
varies with the position in the wire between and 1. 
If the inelastic scattering among the electrons is very 
strong, but negligible between electrons and phonons 
(relevant for short wires), local thermal equilibrium is 
re-established and n £ is a Fermi function with a spatially 
varying electron temperature T a \(r). The resulting tem- 
perature profile can be determined by solving the ther- 
mal diffusion equation with a source term resulting from 
the local energy dissipation. 68 For very long, but narrow 
wires the profile of T e \(r) becomes flat but remains ele- 
vated with respect to the phonon temperature, since the 
electron-phonon coupling becomes very weak under these 
conditions. 



C. Interference of elementary Fermi- and 
Bose-systems 

So far, we have not taken into account quantum in- 
terference between the different modes, or elementary 
Fermi- and Bose-systems. The effect of interference is 
simply that two or more modes are not independent any- 
more, but form new modes, i.e., the coherent superposi- 
tions of the interfering modes. These new elementary 
sub-systems can be charged with energy, entropy and 
particles only as one entity, with a transmission coeffi- 
cient T(e) telling us how constructive or destructive the 
interference is. The simplest example for such a coherent 
superposition are the standing waves, which are formed 
in a quantum well, or a finite quantum wire, by mul- 
tiple reflection at the confining potential walls. In this 
case the new systems cannot support a stationary cur- 
rent. Finite currents can be carried only by propagating 
waves, i.e., scattering 'states', provided that the phases 
of the interfering waves are adjusted such that the trans- 
mission probability is finite. For a completely destructive 
interference again a standing wave is formed, which sup- 
presses the transport of entropy and particles, and the 
transport currents impinging on the interferometer are 
reflected completely. If the source that charges the prop- 
agating modes is at a finite temperature, the local current 
densities of energy, entropy and particles are invariably 
connected. 

The elastic scattering does not produce entropy, but 
results in a coherent branching of the flows of energy, en- 
tropy and particle currentsby the scatterer. This state- 
ment remains valid for arbitrary complex scattering re- 
gions and any number of terminals. The phase coherence 
can be made visible, if the partial waves are brought to 
interference, e.g., by adding a second scatterer or semi- 
transparent mirror, resulting in an electronic or photonic 
interferometer of the Fabry-Perot, 69 Michelson- or Mach- 
Zehnder type. 70 

For these reasons the possibility of phase-tuning of the 
transmission coefficients 7~k, and a corresponding phase- 
dependent thermal conductance 71 not very surprising. 
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The close correspondence between the transport of par- 
ticles and entropy resulting from the present approach 
renders entropy- or 'heat '-interferometer as natural as a 
quasiparticle-interferometer, despite the fact that the en- 
tropy carried by the elementary Fermi- or Bose-systems 
results from an incoherent superposition of states with 
different particle numbers. The entropy stored within 
one elementary Fermi- or Bose-system leaves its coher- 
ent superposition with others unaffected, as the latter 
influences only the transmission coefficients, but not the 
Nk (T, p) and Sk (T, p) responsible for their thermody- 
namic properties. The phase modulation of the entropy- 
or 'heat '-transport by optical interferometers must be 
ubiquitous, but was probably not noted, because the 
heating power transmitted by single optical modes is 
limited by the thermal conductance quantum Lq (see 
Eq. 77), and does not produce large effects at room tem- 
perature. 

Quantum interference also occurs in diffusive systems. 
In this case, the modes of the system cannot be chosen as 
plane waves anymore, but form complex scattering states 
formed by the interference of waves on multiple scatter- 
ing centers. The dominating interference contribution 
to the transport comes from pairs of time reversed dif- 
fusion paths, which return to the starting points. The 
pairs of time reversed paths interfere destructively (at 
zero external field). The interference increases the prob- 
ability of return to a given point (Wreturn = l/V inDt 
for one dimension) by a factor of 2, and thus results in 
a reduction of the conductivity. This effect, know as 
weak localization, can be taken into account as an in- 
terference contribution to the diffusion constant D with 
respect to Eq. 25. Coherent backscattering of light has 
also been observed. An external magnetic field tunes 
the character of the interference between the different 
diffusion paths in a continuous way between constructive 
and destructive (Aharonov-Bohm effect). 54 ' 55 The total 
transmission probability of a set of interfering modes 
strongly depends on the wavelength (i.e., on e) and via 
the Aharonov-Bohm effect on the magnetic field B - lead- 
ing to characteristic conductance fluctuations when ep or 
B are varied. 54,55 

Another type of interference occurs between pairs of 
electrons at similar energies. The Coulomb interaction 
affects the the characteristic energies e of the elementary 
Fermi-systems around the Fermi level, and results in a 
suppression of the DoS near ep. This effect can be taken 
into account as a quantum correction to the particle ca- 
pacity v = dn(T, n)/d[i (Altshuler-Aronov effect). 54,55 
Similar quantum corrections are also expected in the 
thermal conductance, 56 but harder to measure with suf- 
ficiently high precision. 

In conclusion, the concept of elementary Fermi- and 
Bose-systems turns out to be extremely flexible. It can 
be adapted to a wide range of applications in modern 
physics. Once accepted, it provides a more reliable guide 
for our intuition than the classical concept of moving 
particles, as it incorporates the non-classical concepts of 



quantum interference and indistinguishability from the 
start. 



D. Irreversibility and the Loss of Phase Coherence 

The most ingenious side of the Boltzmann equation is 
the fact that the aspect of irreversibility, i.e., the genera- 
tion of entropy is incorporated in the ad-hoc assumption 
of the existence of an inelastic relaxation mechanism, and 
the corresponding characteristic length A; n . As it was 
noted very early, such a mechanism is inconsistent with 
the notion of Hamiltonian dynamics. Any system with a 
discrete energy spectrum is subjected to the recurrence 
objection, i.e., its time evolution must be reversible. The 
recurrence objection is removed by assuming the exis- 
tence of an infinite thermal bath with a continuous spec- 
trum, in which energy and entropy can be dumped with- 
out recurrence. Such an approach is successful for sys- 
tems with a single or a few macroscopic degree of free- 
dom such as quantum bits, coupled to many microscopic 
degrees of freedom. Irreversibility is then generated by 
'tracing out' the bath degrees of freedom. 

Phenomenologically, the generation of entropy can be 
accounted for by damping out the off-diagonal elements 
of the density matrix, which are responsible for the co- 
herent Hamiltonian dynamics. In the limit of long times, 
the decoherence becomes complete, implying a density 
matrix, which is diagonal in the basis of energy eigen- 
states, and with the probabilities {Wj} as eigenvalues. 
Such a diagonal density matrix corresponds to the state 
of maximal entropy under the constraints of Eq. 37, if 
the Wi are given by Eq. 41. The joint dynamics of the 
system and the bath features thermally induced temporal 
fluctuations of physical quantities of the system. These 
fluctuations obey the recently much discussed fluctuation 
theorems. 

A first principles derivation of irreversibility remains 
a severe problem, as the 'first principles' at hand are all 
reversible. In the opinion of the author, it is not clear, 
whether the mathematical operation of 'tracing out the 
bath degrees of freedom' has a correspondence on the 
experimental side. Moreover, there are situations like the 
collisions of heavy ions at very high energy, where vast 
amount of entropy are generated on such short time scales 
(1 fm/c w 10~ 23 s) that there may exist no sufficiently 
strongly coupled bath. 73 



X. CONCLUSIONS 

The purpose of this paper is to provide a coherent 
and self-contained description of the transport of parti- 
cles and entropy both in the macroscopic and the meso- 
scopic regime. To implement this program, the concepts 
of thermodynamics first have to be formulated in a way 
that avoids unnecessary limitations. In order to connect 
the general principles of thermodynamics to the quan- 
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turn physics of matter, the idea of ballistically moving 
particles or plane wave propagation, respectively, has to 
be stripped from all classical elements. I propose to use 
the eigenmodes of the matter field in the language of 2 nd 
quantization as the elementary building blocks of such a 
description. Viewed as thermodynamic systems, called 
here elementary Fermi- and Bose-systems, they can be 
taken as a basis for a unified description of both global 
thermodynamic equilibrium and the ballistic and diffu- 
sive quantum transport. In this description 'classical' 
and quantum transport are the same, while the only de- 
marcation line runs between regimes, where dissipation, 
i.e., the loss of phase coherence, occurs locally, or re- 
motely in macroscopic reservoirs. The notion of elemen- 



tary Fermi- and Bose-systems may prove useful not only 
in solid state physics, including the presently unfolding 
fields of spintronics, caloritronics, and spin caloritronics, 
but also in the description of ultracold atomic and molec- 
ular gases. 74 
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